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Abstract 

We present an unprecedented spectroscopic survey of the Ca ll 
triplet + O I for a sample of 14 luminous (—26 > My > —29), interme¬ 
diate redshift (0.85 ^ z < 1.65) quasars. The ISAAC spectrometer at 
ESC VLT allowed us to cover the Ca ll NIR spectral region redshifted 
into the H and K windows. We describe in detail our data analysis 
which enabled us to detect Ca ll triplet emission in all 14 sources (with 
the possible exception of HE0048-2804) and to retrieve accurate line 
widths and fluxes of the triplet and O l A8446. The new measurements 
show trends consistent with previous lower .2 observations, indicating 
that Ca 11 and optical Ee ll emission are probably closely related. The 
ratio between the Ca ll triplet and the optical Ee ll blend at A4570 A 
is apparently systematically larger in our intermediate redshift sample 
relative to a low- 2 ; control sample. Even if this result needs a larger 
sample for adequate interpretation, higher Ca ll/optical Ee ll should be 
associated with recent episodes of star formation in the intermediate 
redshift quasars and, at least in part, explain an apparent correlation 
of Ca II triplet equivalent width with 2 ; and L. The Ca ll triplet mea¬ 
sures yield signiflcant constraints on the emitting region density and 
ionization parameter, implying Ca II triplet emission from log Uh A 
11 [cm“^] and ionization parameter logU < —1.5. Line width and 
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intensity ratios suggest properties consistent with emission from the 
outer part of a high density broad line region (a line emitting accretion 
disk?). 

Subject headings: quasars: general - quasars: emission lines - lines: 
profiles - black hole physics 


Introduction 


Explaining the origin of Fe emission in quasar spectra is a long-standing 
problem in studies of active galactic nuclei (AGN). Recent observations and models 
fail to reach a consensus on the ionization mechanism even in the case of moderate 
Fe II emitters. Purely collisional models are unable to account for the spectral 


energy distribution of Fe ll emitters (Kuehn et al. 2008). At the same time, the 


ratio UV Fe ll to optical Fe ll emission is more easily explained in the context 


of collisional ionization (Sameshima et al. 2011). Akn 564 Fe ll reverberates in 


response to continuum changes (Shapovalova et al. 2012), providing support for 


the role of photoionization. At the same time Akn 120 Fe ll does not show any 


response (Kuehn et al. 2008). The extreme complexity of the Fe ll ion makes 


theoretical model calculations very difficult and line blending makes estimation of 
Fe II width and strength parameters uncertain. 

Given the difficult interpretation and prediction of the Fe ll spectrum, the 
study of simpler ionic species like Ga ll (emitting the IR triplet A8498, A8542, 
A8662, hereafter GaT) and O l is more straightforward. The Ga ll ion is much 
simpler. The ionization potential of neutral Galcium (?a6.1 eV) leads us to expect 
that Ga ll ions exist where Hydrogen is not fully ionized. Several lines of evidence 
suggest that the GaT and optical Fe ll are produced in the same region. Data 
from Persson (1988) and photoionization calculations ( Joly||1989 ) found that GaT 
is emitted by gas at low temperature (8000 K), high density (>10^^ cm“^) and high 
column density (>10^^ cm“^) similar to optical Fe ll. Matsuoka et al. (2007, 2008) 


computed photoionization models using the O l A8446 and A11287 lines and the 
GaT and found that high density (~ 10^^'^ cm“^) and low-ionization parameter 
iu r\j 10 are needed to reproduce flux ratios, consistently with the physical 
conditions expected for optical Fe ll emission. Density and ionization parameters 


are in agreement with the ones found by Joly (1989). 
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Ferland & Persson (1989) improved the photoionization models including 
physical processes like free-free, H“ bound-free and Compton recoil ionization. 
These authors demonstrated the need for very large column densities {Nc > 
cm“^) to reproduce the Ca ll spectrum. Such large column densities could be pro¬ 


vided by an accretion disk. A similar behavior to the one described by Joly (1989) 


for CaT/H/5 and Fe ll/H/5 was found by Dultzin-Hacyan et ah (1999). These 
authors also suggested that the line emitting region could be associated with the 
outer part of an accretion disk. The width measured on the unblended Fe ll All 126 
line (that belongs to the so-called 1/rm Fe ll lines), is very similar to the ones of 
CaT and O l A8446 and narrower than the one of the Hydrogen lines. In a purely 
kinematical interpretation, this means that low ionization lines like Fe ll, CaT and 


O I A8446 are emitted in the outer zone of the broad line region (Rodriguez-Ardila 


et al.||2002a|b i Matsuoka et al.||2007 , 2008), probably associated with the accretion 
disk, with physical conditions different from the region emitting most of the high 
ionization lines. 


Physical modelling (via a photoionization code like CLOUDY, Ferland et al 


2013) requires input from measurements that are retrieved through model fits of 


the observed spectra. The Ca ll triplet and O l A8446 lines can be easier to model 
than Fe ll but they are not as strong as Fe ll optical features which, however, 
are heavily blended. They are located in a wavelength range where the putative 
accretion disk continuum and the high energy tail of hot dust emission form a sort 
of valley in the spectral energy distribution. The old stellar population associated 
with the quasar host galaxy peaks at ~ Ifim (Landt et al. 2011). In many low- 


luminosity quasars CaT and O l A8446 are confused with the underlying absorption 
of the host galaxy. 

The near-infrared region (NIR) has not been easy to observe due to the lack 
of high quality IR spectrometers. Only in relatively recent times the development 
of new instruments has made it possible to obtain high S/N spectra for sources 
with redshift larger than z > 0.1. It is not surprising that the Ca ll and O l 
features remain relatively unstudied even in sources for which a wealth of optical 
data exist. 

The present work extends the study of the Ca ll triplet in emission to include 
14 high luminosity and intermediate redshift quasars ()p|). Whenever possible they 


are analyzed within the 4D eigenvector 1 context (4DE1; Boroson & Green 1992 


Sulentic et al. 2000a|b i Marziani et al. 2001, 2003a|b ; [Sulentic et al. 2007), which 
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serves as a spectroscopic unifier/discriminator of the emission lines properties for 
type 1 AGN. The parameters defining the four dimensions of the El space involve: 
(1) FWHM(H/3 b( 7), (2) the ratio between the equivalent width (W) of the optical 
Fe II A 4570 blend and H/5, i?Feii= W(Fe ll A4570)/W(H/3) /(Fe ll A4570)//(H/5), 
(3) the soft X-ray photon index (Fsoft), and (4) the centroid line shift of high- 
ionization C IV A1549, with the strongest correlations involving parameters 2, 3 
and 4 (Sulentic et al. 2007j ). In the 4DE1 framework the broad line AGN can 
be divided into two populations, A and B (Sulentic et al. 2002). Gonsidering 
the broad component of the hb line, populations A and B can be separated at 
FWHM(H/5sc) = 4000 km s“h Spectra of sources above and below this limit 
look substantially different. It is possible to introduce a finer subdivision that will 
not be used in this paper because of its small sample size: we will distinguish 
sources as Pop. A or B only. Quasars can be divided into bins of AFWHM(H/5) 
and Ai?Feii- The bins of Pop. A (Al, A2, A3, A4), are defined in terms of 
increasing Rpeii, while Pop. B bins (Bl, Bl"'', Bl'*'''') are defined in terms of 
increasing FWHM(H/3). Population A shows: (1) a scarcity of RL sources, (2) 
strong/moderate Fe ll emission, (3) a soft X-ray excess, (4) high-ionization broad 
lines (HIL) with blueshift/asymmetry and (5) low-ionization broad line profiles 
(LIL) best described by Lorentz fits. Population B: (1) includes the large majority 
of RL sources, (2) shows weak/moderate Fe ll emission, includes sources (3) with 
less prominent or no soft X-ray excess ( [Sulentic et al.| 2007) and (4) with HIL 
blueshift/asymmetry or no blue shifts at all. Last, (5) Pop. B shows LIL Balmer 
lines best fit with double Gaussian models. The physical drivers that change along 
the 4DE1 sequence have been identified: number density appears to increase from 
Pop. B to A, and black hole mass increases with large scatter from A to B. The 
principal driver of source occupation in the 4DE1 space is probably Eddington 
ratio which increases from B to A (Marziani et al. 2001; Boroson 2002). The 
4DE1 parameter space is, to a first approximation, not strongly dependent on 
luminosity. The same basic distinction between Pop. A and B is recognizable also 
for very luminous sources (Marziani et al. 2009). 


This paper presents new observations and data reduction of spectra of the 
Ca II IR triplet, as a representative LIL, as well as of the O l A8446 line, in very 
luminous quasars of the Hamburg-ESO survey ^and Q. We discuss in some 
detail the identified line and continuum components for a proper data analysis (^. 
Basic results obtained from the measurements are briefly reported in In ^we 
compare observational data and photoionization predictions, the role of CaT and 
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O I A8446 within the 4DE1 context, possible implications for star formation, and a 
preliminary geometry for the broad line region (BLR). Conclusions are presented 
in ^ 


2. Sample selection 


Our original sample of Hamburg-ESO (HE) sources was defined with the aim 
of testing/extending the 4DE1 correlation at relatively high- 2 ; and at high lumi¬ 
nosity ( [Marziani et al.|[2009| . The HE original sample of H/3 ISAAC observations 
was therefore flux limited and unbiased in terms of Fe ll strength. Our Ca ll 
subsample was extracted from the 52 sources ISAAC sample of [Marziani et ah 
(2009). The only selection criterion was that IR atmospheric absorptions were 


not affecting the Ca ll and O l lines. All targets in the present investigation have 
therefore optical spectra around H/3 already obtained with the same instrument. 

The Ca ll -|- O l sample contains 15 high luminosity Hamburg-ESO (HE) 
quasars, with Mb < -26, in the redshift range 0.846 < 2 : < 1.638. They were ob¬ 
served at the Very Large Telescope (VLT-UT3) equipped with the Infrared Spec¬ 
trometer And Array Camera (ISAAC) that yielded a spectral resolution Rs ~1000. 
Sample characteristics and the log of the observations are reported in the Table 
The first column lists the coordinate name from the Hamburg ESO survey for each 
object of our sample. In columns 2, 3, 4 and 5 we report redshift, apparent and 
absolute magnitude, and Kellerman’s for radio-loudness parameter, as previously 
given by Sulentic et al. (2004) and Marziani et al. (2009). Column 6 lists the 


population in the Eigenvector 1 sequence (Sulentic et al.|2002). One of the sources 


(HE2259-5524) was excluded from the analysis because of poor S/N ratio. 

The median absolute magnitude of the selected 14 HE sources IS -28, which 
makes them among the most luminous quasars known. The sample presented in 
this paper is small; however, Ca ll and O l observations of quasars have been 
until now obtained for only 66 sources (Persson||1988; Matsuoka et al. 2005, 2007 


2008; Landt et al.||2008 ), of which only one at 2 : > 1. Therefore the present sample 
represents the first systematic study of the Ca ll and O l lines in quasars at redshift 

> 1 . 


In the following we will consider as a reference low- 2 :, low-luminosity “con¬ 
trol” sample the sources of Persson (1988) that were observed with very high S/N 
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and spectral resolution similar to the one achieved in the present paper. This 
previous Ca ll sample was sources with strong Fe ll emission at low redshift, as are 
other, more recent ones (Matsuoka et al. 2007). We stress again that our sample 
is not selected on the basis of Fe ll intensity. The selected ISAAC sources do not 
show very strong optical Fe ll emission: 10 of them are Pop. B and 4 Pop. A (Ta¬ 
ble [^. The Pop. A sources are of spectral type Al following Sulentic et al. (2002); 
for all sources we therefore have i?Feii ^0.5. This allows us to explore a domain 


in Fe II emission and in the 4DE1 parameter space where the study of Ca ll could 
give new insights on the origin of low-ionization lines for the majority of quasars 
with moderate Fe ll emission. At the same time, the preferential selection of high 
i^Feii sources at low-^: should be taken into account when possible luminosity and 
effects are studied. 


3. Observations 


Observations were aimed at obtaining spectra of the O l A8446 and CaT 
lines for the sample described in the previous section. Considering the rather low 
efficiency of IR spectrometers available until a few years ago, and the relatively 
low equivalent width of the features we wanted to observe (~ 10 — 20 A), the 
collecting area of a large telescope was needed to obtain sufficient signal with a 
moderate dispersion. 


The VLT includes four telescopes of 8.2 m diameter located on top of Cerro 
Paranal, and instruments are mounted at the foci of the four telescopes. ISAAC 
wa^one of them, and has been for many years one of the few instruments available 
worldwide for moderate resolution IR spectroscopy of faint sources like high-^; 
quasars. ISAAC was able to obtain images and spectra in the wavelength range 
1-5 /rm ( Moorwood et ah] 1998). It was equipped with gratings for high (i?s ~ 
3000) and low spectral resolution (Rs ~ 1000). 


Our spectra were collected on the Rockwell CCD detector (pixel size 18.5 /rm 
yielding a pixel scale of 0.148 ’’/pixel) at low resolution but with a slit width of 
0.6”, that ensured R^ ~ 1000, as measured on sky and arc lamp lines. The second 
part of Table provides date of observation, photometric band of the covered 
ranges, detector integration time (DIT), and number of sequential exposures with 


ISAAC has been decommissioned in 2013. 
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fixed DIT in the spectral band covering the CaT + O l blend. Seeing values, 
reported in the penultimate column of Table [T| are averages over 10 - 15 measures 
at Paranal obtained during each quasar exposure. The S/N was measured on the 
continuum at 8300 to 8400 A. The instrumental spectral resolution at FWHM 
measured on sky lines was ~ 330 km s“h 

4. Data reduction 

We used the IRAF software to reduce our spectra. Frames were obtained 
alternating exposures of two or three DITs in two different positions (A and B) 
along the slit, following a sequence like A-B-B-A. The frames in each position 
were averaged and then subtracted one from the other. This procedure cancelled 
the bias and dark current present on each frame, assuming that they remained 
constant during the sequence, and also subtracted most of the sky background, 
facilitating the spectrum extraction. After the subtraction, frames were divided 
by the corresponding flat fields. We extracted the ID spectra using the iRAF 
task APSUM, where a low-order polynomial function was fitted and subtracted to 
eliminate any residual background. ID spectra were wavelength calibrated using 
an Ar-Xe lamp. The calibration was achieved using a 3rd order polynomial, with 
a typical rms of 0.3 A in the H band and 0.4 A in the K band. This calibration, if 
applied to the sky spectrum yielded a slight displacement that was corrected using 
measures of OH sky line wavelengths. We applied the offset-corrected wavelength 
calibration to the quasar spectra. Finally ID quasar spectra were combined to get 
one median spectrum. The spectra of the standard stars (observed right before or 
after the quasar spectrum) were extracted and wavelength calibrated in the same 
way. 

We corrected for the effect of telluric lines using the TELLURIC IRAF task. 
This procedure involved a change in the quasar continuum shape since the task 
TELLURIC carried out a division between the quasar and the standard star. To 
return to the original quasar continuum after the elimination of the telluric bands, 
we divided the quasar spectrum by the standard star continuum (i.e., by the 
standard star divided by a synthetic atmospheric transmission spectrum). 

Fourteen standard stars were spectral type B and one G. The B stars belonged 
to different luminosity classes (supergiants, giants and main sequence), while the 
star G is a main sequence star. Comparing the NIR stellar spectral library from 


Ranade et al. 

(2004 

) and 

Ranade et al. 

(2007 

) with theoretical libraries in the H 

and K bands ( 

Pickles|1998 

), we found that the theoretical spectra could reproduce 


the observed spectra of main sequence stars, but not those of giants and super¬ 
giant stars. For the giants and supergiants stars we considered a spectral energy 
distribution provided by black-body emission at the temperature tabulated for the 


star (Ammons et ah 2006). Black-body emission fits very well the spectra after 


erasing the stellar absorptions lines with the tool xydraw of the SPLOT iraf task 
(see Figure [^. In order to render our approach uniform we fitted all the stars 
using a black-body, because there is practically no difference between the simple 
black-body fit and the continuum spectra of the stars. Finally, the standard star 
spectra were normalized to H and K 2MASS magnitudes. 

The sensitivity function was obtained dividing the standard star spectra by 
the black-body model. The final relative flux calibration was achieved dividing 
the spectrum by the sensitivity function. Spectra were corrected for Galactic 


extinction following Schlegel et ah (1998) but no internal extinction correction 


was applied. Quasar spectra were normalized to 2MASS magnitudes as done for 
the standard stars in order to minimize the effect in the continuum flux change 
possibly induced by intrinsic variability (see next section). Redshift correction 
was performed using the redshift values reported in [Sulentic et ah (2004) and 
Marziani et al. ( 2009|| ). These redshifts were measured on H/5jvC) Hyjvc and 03 
with an uncertainty usually cz < 150 km s“^ in the rest frame. Figure shows the 
rest-frame flux and wavelength calibrated spectra before continuum subtraction. 


4.1. Variability 


The H/3 spectra were taken between 2001 and 2005, while the near-infrared 
spectra were collected in 2010. There is a difference of 5 - 10 years between the 
H/3 and NIR observations. In order to estimate a possible broad line flux variation 
in response to continuum changes, we used the BLR radius vs. luminosity relation 


of Bentz et al. (2013). Comparing the expected response time with the difference 


in date between the optical and IR observations, we found that eight objects could 
have varied. Extracting the light curves from the Catalina Real-Time Transient 
Survey (CRTS; Drake et al.||2009 ), we could identify eight cases of actual variation. 
Five of the eight objects show a systematic variation <0.1 mag in a period of 6 
- 8 years when they were monitored by the CRTS. The remaining three objects, 
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HE0048-2804, HE0203-4627 and HE2147-3212, show a “high frequency” variation 
which is apparent even from a visual inspection of the light curve. The amplitude 
is modest, 0.2 magnitudes, and does not hint at any systematic variation over 
longer periods. We normalized the Ca ll spectra to H and K 2MASS magnitudes, 
and the H/3 spectra to J and H magnitudes, in order to obtain an accurate flux 
scale as well as to remove any variability effect on line and continuum fluxes. This 
normalization was not possible for two objects, HE0048-2804 and HE2340-4443, 
since 2MASS does not have information on the I or Z photometric bands that 
cover redshifted H/5 (see Appendix [A| for a brief discussion). 


5. Measurements 


The flts and data analysis were carried out within the Eigenvector 1 context 
(Boroson & Green 1992; Sulentic et ah 2000a|b i Marziani et ah 2001, 2003a|b 


Sulentic et al. 2007). We performed multicomponent flts including all the rele¬ 


vant components that we presume are contributing to the spectral ranges under 
investigation. The objects analyzed in this paper belong to a sample of intermedi¬ 
ate redshift quasars for which H/3 is available from previous ISAAC observations. 
The H/3 line profile parameters can guide us in the interpretation of the NIR re¬ 
gion. The optical spectra were previously analyzed by Sulentic et al. (2004) and 


Marziani et al. (2009). The analysis of this paper adds an H/3 multicomponent 


SPECFIT modelling that was not carried out for individual sources or, if carried 
out, was not presented in previous papers. The optical and the NIR spectra were 
modeled with the task SPECFIT of IRAF, a routine that carries out simultaneous 
flts of several line and continuum components to ID spectra (Kriss 1994). Each 
component is described by a set of parameters that are varied according to the 
initial guess between a maximum and a minimum value. The best fit is achieved 
via minimization. The multicomponents fits are shown in the Figure The 
left panel shows the optical spectra and the right panel the NIR spectra. 


Few parameters were actually free to vary in the NIR spectral range. The H/5 
profile was taken as a proxy for the Pa9 A9229 and O l A8446, and the peak shifts 
of these two lines were set as measured on H/3. The flux ratio between the broad 
component and the very broad component of Pa9 is assumed to be the same of 
H/3. The width and shift of 03 set the width and shift of [S ill] A9531, so that 
only the intensity of the latter line is left free to vary. The centroid of CaT was 
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set in the quasar rest-frame. The only remaining free parameters were therefore: 
the specific fiux scale and the index of the power-law continuum (2), the intensity 
scaling of the Fe ll template (2), the intensity and width of O l A8446 and CaT (4), 
the intensity scaling of the high order Paschen lines (1) and the intensity and tem¬ 
perature of Paschen continuum (2), giving a total number of 11 free parameters. 
However, not all these parameters refer to features blended with O l A8446 and 
CaT: the Paschen continuum properties do not affect the blend O l -|- CaT; the 
intensity scalings of the high order Paschen lines and of the Fe ll template are also 
set by features outside of the blend. This makes the fits more robust, and ensures 
that O I A8446 -|- CaT intensity and width are the parameters mostly affecting the 
xl in the wavelength range they cover. The following sections describe each one 
of the components used for the fits and the error estimates. 


5.1. Continuum 


The main continuum emission present from the UV to the NIR spectral re¬ 
gion of quasars is due to an accretion disk ( Malkan &: Sargen^|1982 ; Malkan||1983 ). 
Accretion disk emission is thermal (a “stretched black-body”) but can be modeled 
locally with a power law with a slope of —2.4. Optical and NIR spectra should 
share in part the same underlying continuum i.e., the low-energy tail of the ac¬ 
cretion disk emission. However, we could not consistently fit the continuum with 
a single slope from the optical to the NIR. In the majority of the sources the NIR 
spectrum is ~ 30 - 80% above the NIR extrapolation of the optical continuum. 
Therefore, a local power-law continuum was fitted for each quasar assuming a free 
power-law index. The optical continuum was defined at 4750 and 5100 A. The 
IR continuum was set at 8100, 8800 and/or 9400 A, depending on the wavelength 
range covered by each NIR spectrum. 

A possible origin of the disagreement between optical and NIR continuum 
is a contribution from two additional NIR components (Landt et al. |2011 ). The 
first is a diffuse continuum (due to thermal emission and scattering) that may 
flatten the NIR spectrum ( [Korista fc Goad 2001). The second component is due 
the finite albedo of the same dusty clouds emitting the thermal far-IR continuum 
(Landt et al. 2011). The small range covered by our spectra makes it impossible 


to measure their contributions. The continuum can be also affected by the host 
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galaxy (this seems to be case for only one quasar in our sample, see Section 5.8). 


5.2. Fe II template 


Over the years, there have been several efforts to model the ultraviolet and 


optical Fe ll contribution (Sigut & Pradhan 2003; Sigut et al. 2004; Bruhweiler 


fc Verner] 2008). In the NIR, a semi-empirical Fe ll template based on the I Zw 


1 spectrum was produced by Garcia-Rissmann et al. (2012) who kindly made it 
available to us. A theoretical template based on photoionization was also computed 
by the same authors (Table 3 of Garcia-Rissmann et al.|2012 ). A comparison of the 
two templates is shown in the upper left panel of Figure]^ The agreement between 
the two templates is fair; both indicate significant Fe ll emission at around 8400 
A, 9200 A, and l/rm. These features are expected to be produced through Lya 
pumping of high energy levels (j^i 11 eV) in the Fe"'' ion, although Lya pumping 
may not be the only production mechanism, and additional mechanisms are still 


debated (Sigut & Pradhan 

1998; 

Rudy et al. 

2000; 

Rodriguez-Ardila et al. 

2002a; 

Garcia-Rissmann et al. 

2012 

). Lya fluorescence makes the difference between 


optical and NIR emission. At present we do not have evidence that this process 
could induce a difference between the NIR Fe ll spectrum of Pop. A and B sources. 
We considered Fe ll measures of Table 6 of Landt et al. (2008) for Pop. A (13) and 
B (4) sources. The median Fe ll relative intensities measured in the range 8000 - 
9810 A are the almost identical. 


The upper right side of Figure]^ shows an example of fits to the same spectrum 
(HE1349-I-0007) using different templates. The lower panels of Figure]^ show the 
HE1349-I-0007 NIR spectral range after background subtraction, to illustrate the 
effect of the Fe ll template on O l A8446 and GaT measures. 

In our spectra with high S/N we observe a rather symmetrical bump due to 
Pa9 -|- Fe ll at ~ 9200 A. The emission on the blue side of this bump cannot be 
reproduced with the semi-empirical template. The width of Pa9 is limited by the 
assumption that it should be consistent with the one measured for the Balmer lines 
in the optical spectra. For this reason we could not model Pa9 with a profile broad 
enough to fill the observed excess emission. In addition, higher order Paschen lines 
constrain the Pa9 intensity. Therefore, the excess of emission with respect to the 
semi-empirical template has to be ascribed to Fe ll. Given the problem with the 
semi-empirical template, the Fe ll emission of our spectra was fitted with both the 
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semi-empirical and the theoretical template. 

Between 8600 and 8800 A the semi-empirical template presents a contribution 
that is not modeled by the theoretical one. This feature affects the flux and 
equivalent width of the line Ca ll A 8662 A. Because the three lines have the same 
intensity, if the semi-empirical template is used, the equivalent width and flux of 
the CaT and O l A8446 decrease. In our sample the O l A8446 and CaT flux and 
equivalent widths is affected by less than 20%, save for: two sources where the 
difference reaches 20 - 35% for O l A8446 (HE2340-4443 and HE0203-4627), 
and four sources where the CaT differs by 20 - 35% (HE0203-4627, HE2352-4010, 
HE0248-3628, HE2340-4443). 

We then choose the template that gave the minimum xl value. We found 
that the theoretical template suits better in the majority of the cases. In only 
three cases (HE2147-3212, HE2349-3800 and HE2352-4010) the semi-empirical 
template was used; in the remaining sources the theoretical template was applied. 

In the IR, the strongest, most isolated Ee ll emission features are around 1 
/rm, but this wavelength range is not covered by our spectra. In the wavelength 
range covered by our spectra we have no truly isolated Fe ll features. However, we 
could still accurately measure the Fe ll contribution in the range of 9000 - 9400 
A from the scaled and broadened template, and consider it as an estimate of the 
NIR Fe II emission. In the majority of our spectra, we were able to neatly fit the 
Fe II and Paschen 9 line contribution. 


For the Fe ll optical contribution in the H/5 region we used the template 
previously employed by Marziani et al. (2009), mainly based on a high S/N optical 
spectrum of I Zw 1. Iron emission gives rise to a fairly isolated bump at ~ 4450 
- 4680 A. Optical Fe ll emission was defined by the total emission in this range 


(conventionally indicated with Fe llA4570), a standard practice (Boroson & Green 

T^. 


5.3. Broad component 


Following the 4DE1 approach (Sulentic et al. 2000a|bi Marziani et al. 2001 


2003a||b ; Sulentic et al.|2007 ), we assume that low-ionization lines in Pop. A and B 


sources have profiles that are best fitted by Lorentzian (Pop. A) or double Gaussian 
(Pop. B). All the broad components were modeled with this assumption. We note 
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that the Lorentzian approach is roughly equivalent to fitting two Gaussians with 
the same shifts (usually < 200 km s“^ for Pop. A sources, Sulentic et al.|20l'2 ), one 
broader and one narrower. However, this is strictly true if the S/N is low: at very 


high S/N (not the case of our data) a Lorentz function yields a lower xl (Marziani 
et al.|[^03b ). Since the best data suggest that Pop. A H/5 profiles are Lorentzian 
(Veron-Cetty et al. 2001; Marziani et al. 2003b| Zamfir et al. 2010; Shapovalova 
et al. 2012), we fit Lorentzian functions to the broad line profiles of all Pop. A 


sources in our sample. With a double Gaussian fit, the difference between Pop. A 
and B would be that the broader Gaussian is usually redshifted in Pop. B. 

The H/3 line profile was taken as a reference for modelling the GaT and O l 
A8446 blend, since its profile is well defined and well understood. GaT was modeled 
with three BGs of equal intensity and the same FWHM since it is believed that 
the Ga ll triplet components are optically thick (Persson||1988; Ferland & Persson 


1989). As Ga ll lines are completely blended, it is difficult to measure a shift, 


so that the Ga ll lines were held fixed at rest-frame wavelength. The maximum 
shifts possible for O l A8446 were the ones measured on H/3, due to the relation 


between the two lines (see also section 7.1). FWHM and shifts of H/3 and high 


order Paschen lines were assumed to be identical within the uncertainties. 


5.4. Very broad component 


A VBC was fitted to H/3 and the most intense high order Paschen lines (Pa9, 
Pae A9546 and Pa5 A10049) for all Pop. B sources. The shift and FWHM were 
assumed equal to those of H/5 within the uncertainties. A VBG was not fit to Ga 
II triplet, since the VBG is believed to be associated with a very high ionization 


region (the very broad line region, VBLR e. g. Snedden & Gaskell 2007 and 


Marziani et al. 2010) that is expected to emit negligible Fe ll and Ga ll. Since 


there is no clear spectroscopic evidence of the presence of a VBG for O l, only a 
BC was first considered also for this line (although the possibility of an O l VBG 


is debatable and discussed in section § 7.2). 


5.5. Narrow components and narrow lines 

We considered the H/5 narrow component (H/3 nc) and 03. In the NIR, we were 
able to detect a hint of the NC for O l only in HE1349-I-0007, a source in which 
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H/5nc is especially strong. Also, in all the spectra covering the region around ~ 
9500 A, we could detect [S ill] A9531. 


5.6. High order Paschen lines 


Modelling the NIR spectral range with only the strongest high order Paschen 
lines visible in our spectra (Pa9 A9229, Pae A9546 and Pa(5 A10049) results in an 
emission deficit between 8700 and 9100 A. Therefore we decided to include even 
higher order Paschen lines: from Pa9 A9229 up to Pa24 A8334. Since Pa9 is clearly 
seen in the majority of our spectra, it was taken as a reference to scale the intensity 
of the other Paschen lines. The intensity ratios were computed from case B values 
from CLOUDY simulations (version 08.00; Ferland et al.||1998 , 2013) which have an 
exponential decay with increasing order number. This set of lines forms a pseudo¬ 


continuum that cannot be neglected in some cases. Previously, Persson (1988) 
found that Paschen lines contribute ~ 12% of the flux of Ca ll A8662 in Mrk 42, so 
that they were ignored. Comparing the flux of Pal3 A8665 and the Ca ll A8662, we 
found that Paschen lines do contribute by a fraction gtrsim 30% in ten sources; 
in two sources the contribution is ~ 10% (HE0203-4627, HE2147-3212). In two 
cases (HE0048-2804 and HE0058-3231) the high order Paschen lines dominate the 
fit. 


5.7. Paschen continuum 


We could detect a hint of Paschen continuum (PaC) in HE0048-2804 and 
HE2202-2557, so we included an optically thin model of the Paschen continuum 
at 8204 A. The inclusion has little effect on the O l A8446 and CaT measures, 
since the Paschen series head is separated in wavelength from the O l -|- CaT 
blend. A CLOUDY simulation with log U = —2.5 and logun = 12 provided us with 


predictions on the relative intensity of Pa9 and of the integrated PaC (c.f. Oster- 
brock &: Ferland||2006 ). We then estimated the continuum specific intensity at the 
Paschen edge, and assumed an exponential decrease toward shorter wavelengths 
appropriate for an optically thin case. Our estimates appear to be in excess with 
respect to the observations since the best fits require a PaC smaller than that 
predicted. There are three main possible explanations: (1) CLOUDY computations 
over-predict the recombination continuum (a distinctly possible case since the as- 
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sumed conditions are unlikely to be found in optically thin gas), ( 2 ) the quasar 
continuum level is not correct, and the actual continuum is lower, or (3) PaC is 
not optically thin. We cannot test these options on our data since the spectra 
cover only a small wavelength range. To properly dehne PaC we need to cover a 
spectral range spanning from the optical to the NIR. We will try to address the 
issue of the PaC intensity in the future. 


5.8. Contribution of the host galaxy 


A new stellar population synthesis model (Chen et ah in preparation), based 
on the code of Bressan et ah (1998) with updated stellar evolutionary tracks (|Bres 


san et al. 2012 ) and with stellar atmospheres was used to account for the possible 


contribution of the quasar host galaxy. These computations are suitable for the 


analysis of stellar absorption lines in the optical and NIR spectral regions (Sanchez- 


Blazquez et al. 2006b|[cpii Rayner et al. 2009). From the lookback time at the 


quasar redshift we estimated an upper limit to the host age. Then we computed 
the black hole mass, and estimated the host mass assuming the Mbuige/AfRU ratio 


(Magorrian et al. 1998; Merloni et al. 2010) appropriate for the quasar redshift. 


The 0.9 pm specific flux was then derived from the stellar population synthesis 
models that are provided as a function of mass. We found that the underlying 
stellar absorption of the host galaxy is significant only in HE2202-2557, with a 
luminosity contribution of ~ 50%, while the rest of the sample is affected only by 
< 10%. In the HE2202-2557 case, we assumed a spheroid mass ~ 1.1 • IO^^Mq, 
an age of 2.4 Gyr and a metallicity of 2 Z 0 . These properties are consistent with 
massive ellipticals expected to host very luminous quasars at intermediate redshift 


(Trager et al. 

to 

0 

0 

0 

Kukula et al. 

2001 ; 

Sanchez et al. 

2004 

Kotilainen et al. 

2007; 

Spolaor et al. 

2009; 

Floyd et al. 

2013) 


In order to test the behavior of several stellar templates, we performed several 
models taking different spheroid masses ( 10 ^°, 10 ^^, 10 ^^, 10 ^^ Mq, but keeping 
metallicity equal to 2 Z 0 and host age 2.4 Gyr). Model spectra with different 
spheroid masses give a similar total spectrum, but spectra differ in the CaT range 
where the stellar absorptions have different depths and widths, as shown in Figure 
The CaT intensity in absorption decreases with increasing galaxy mass and/or 


velocity dispersion (Cenarro et al. 2003; van Dokkum & Conroy 2012 ) while its 


width increases. These changes affect the emission lines measurements. The flux 










































































16 


and equivalent width of CaT in emission and of O l A8446 increase with galaxy 
mass. The effect induced by a 3 order of magnitude mass change is, in general, < 
50%. In the HE2202-2557 case, Figurej^shows that host masses lower than 1 • 10^^ 
Mq would yield significant residuals and a worse while for 1 • the xt 

will not change significantly. We assume a bulge mass 1.1 • 10^^ Mq, consistent 
with the value expected from Mbuige/Tfeu ratio at intermediate redshift following 


Merloni et al. (2010). 


6. Results 

Figure shows the multicomponent fits after continuum subtraction for the 
H/5 and O l A8446 + CaT spectral regions. Results of multicomponent fits are 
provided in Tables and Table [^reports H/5 and Pa9 A9229 measurements. The 
first column lists the name of the object and the second column is the continuum 
flux measured at the line’s rest-frame. The next columns list equivalent width 
(W), flux (F) and full-width at half maximum (FWHM) for the broad component 
(BC), narrow component (NC) and very broad component (VBC). The last column 
reports the FWHM of the full profile i. e., the sum of BC and VBC. Table [^reports 
the O I A8446 and CaT A8498, A8542 and A8662 measurements. As for Table 
the second column lists continuum flux at the O l A8446 wavelength and the other 
columns report W, flux and FWHM for the BCs. The NC was detected for O l 
only in the case of HE1349+0007. Error estimates were obtained as described in 
Appendix |B| The equivalent width and flux for optical and NIR Fe ll are reported 
in Table The last column reports which NIR template (i.e., theoretical or 
semi-empirical) was used in the fits. Shifts are, at a first guess, assumed to be 
consistent with those of H/3, and are found to be in general close to 0 (i.e., consistent 
with the quasar rest frame). They are therefore not reported. Only in one case, 
HE1409-I-0101, the H/5bc shows a significant blueshift, -700 km s“^, and a similar 
shift was applied to the O l peak to make it consistent with H/3. The Ca ll triplet 
was assumed to be always at rest frame. The VBC of Hydrogen lines was fitted 
with a broad Gaussian that has a large shift to the red, ~ 1000 - 3000 km s“\ 
assumed to be the same for optical and NIR Hydrogen lines. 
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6.1. Ca II triplet detection 


The Ca ll triplet has been observed in sources with strong Fe ll emission, as 


in the majority of the Persson (1988) sample. CaT could be observed in many 


Fe II strong sources because they are Narrow Line Seyfert 1, and emission lines 
are therefore relatively narrow (FWHM H/3 < 2000 km s“^) and sharply peaked. 
In our sources the CaT lines are often broad and totally blended with O l A8446 
and high order Paschen lines making it necessary to confirm that CaT is really 
detected. The A8498 and A8542 fines do not provide conclusive evidence, because 
they are heavily blended among themselves and with O l A8446. On the converse, 
Ca II A8662 is in a region relatively free from strong contaminating lines, and 
in some cases this line could be detected by a simple visual inspection of the 
spectra (HE0248-3628 is an example). In two cases of sources with weak CaT, 
like HE0048-2804 and HE2202-2557, a visual inspection is inconclusive. To test 
that detections are real, we performed a fit with no CaT emission. In HE0048-2804 
we could observe a slightly positive residual around 8600 A; the CaT emission is 
very weak, if any (Eigurej^ left). The positive residual is strong in HE2202-2557 
and without any doubt CaT emission is present (Eigure right). We conclude 
that Ca II triplet is detected in all of our sources with the possible exception of 
HE0048-2804, even if detection of CaT emission was not a priori expected on the 
basis of Ee ll - Ca ll correlation ( Persson||1988| |Joly||1989 i Eerland &: Persson||1989 


Dultzin-Hacyan et al. 1999). We also note that W(CaT) can be as large as ~ 50 

A. 


6.2. Relations between Ca ii triplet and O i A8446 


The equivalent width distributions of CaT and O l A8446 for Persson (1988) 


and our sample are shown in the upper panels of Figure The spread of O l 
A8446 and CaT values is similar in the two samples. In our sample W(CaT) is 
distributed over a wide range, whereas the W(0 l) distribution is peaked around 
~ 15 - 20 A. The Persson sample similarly shows a narrower distribution of W(0 l 
A8446). The difference in the equivalent width distributions is reflected in the wide 
range of the CaT/0 l A8446 ratio: from almost 0 to CaT significantly stronger 
than O I A8446. 


If we add the 11 upper limits in the sample of Persson (1988) to the distri¬ 


bution, the CaT seems to be significantly more prominent in our sample. Two 
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sample tests that include censored data indicate a confidence level between 0.92 
(Gehan’s Generalized Wilcoxon test) and 0.98 (Logrank test). The O l A8446 dis¬ 
tributions (no upper limit; lower left) are instead consistent. |Persson| ( |1988[ ) fails 
to detect CaT emission if GaT/0 l A8446 < 0.3, with W(0 l A8446) ~ 10 - 20 
A. These sources all show i^Peii ^ 0.5, and are of spectral type B1 or Al. This 
finding leaves open the possibility that there could be a dependence on L and/or 
( ^6.5[ ) specially for GaT. 


6.3. Line widths of H/5, O l A8446 and Ca ll triplet 


Persson’s and our data are complementary because our sample is mainly made 
of Pop. B sources (Table column 6). Twelve of the Persson (1988) sources are 
Pop. A for which FWHM(H/3) < 4000 km s“^ and 2 are Pop. B. (ll Zw 136 and 
Mrk 376). If we compare the FWHM of H/3 (measured on the full profile, without 
separating VBC from BC; top panel of Figure we see an overall consistency 
of FWHM(H/3) with FWHM(CaT) and FWHM(0 i A8446) for Pop. A sources, 
while for Pop.B sources FWHM(H/3) is significantly larger than the FWHM of the 
NIR lines. This is the effect of the VBG that broadens the full profile of Pop. B 


by 20 - 30% (Marziani et al. 2013a). 


If we restrict the attention to the broad component of H/3, there is an ob¬ 
vious correlation between FWHM(H/3 bc) ^ FWHM(CaT) ~ FWHM(0 l A8446) 
(middle panel of Figure]^, which confirms and extends the trend found by Pers-| 
(1988[). In the middle panel of Figure]^ Pop. A sources are located around 


son 


the equality line, while Pop. B sources are located below. The Pop. B trend is 
significant at a 2-a confidence level accord into to a Wilcoxon signed rank test 


(Wilcoxon 1945), and may indicate that there is a systematic difference between 


the widths of H/3 BC only and O l A8446 or CaT in Pop. B sources. A discrep¬ 
ancy could also be due to different quality of data and/or to different measurement 
techniques, since Pop. A sources are mainly from the Persson sample and Pop. 
B sources are from our sample and the H/5 profiles need a BC/VBC decomposi¬ 
tion. If broadening is predominantly due to virial motions, and if FWHM(H/5s(7) 
is really systematically larger than FWHM(CaT) and FWHM(0 l A8446), the H/3 
emission associated with the BC only may still preferentially come from regions 
closer to the central black hole than CaT and O l A8446. 


It is also interesting to compare the FWHM of O l A8446 and CaT (bottom 
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panel of Figure]^. At least to a first approximation, O l A8446 and CaT show 
consistent FWHM values. The Persson data suggest that CaT is slightly system¬ 
atically broader than O l A8446. This effect is not seen in our data that have 
however much larger uncertainties. A small CaT and O l A8446 FWHM difference 
may be explained in the context of the different emissivity of the two lines as a 
function of ionization parameter and density. If the two fines are not emitted in 
a coexistence region (Rodriguez-Ardila et al. 2002a; Matsuoka et al. 2005, 2007), 
they can be emitted in two regions with different physical conditions, but kinemat¬ 
ically coupled. Considerations on ionization parameter and density may support 
this possibility (discussed in §7.1). 


6.4. Relation to Fe ii 


Photoionization models performed by Joly ( 1989[ ) have shown that the relation 
between the ratios CaT/H/3 and Fe llA4570/H/5 provides evidence of a common 
origin for CaT and Fe llA4570: CaT/H/5 increases at high density and low tem¬ 
perature as does Fe llA4570/H/3 (Joly |19^ . Our sample follows the same trend 
found by Joly (Figure]^, with a Spearman correlation coefficient p ~ 0.73, yield¬ 
ing a negligible probability {P < 10“^) that the correlation arises from statistical 
fluctuations. The correlation is consistent with the idea that Fe llA4570 and CaT 
are both emitted in a similar region, since the two fines are expected to be emitted 
under very similar physical conditions. If not, a hidden variable (metallicity or 
even orientation: Dultzin-Hacyan, Taniguchi & Uranga 1999) could give rise to 
the observed correlation. 


A least-squares fit including censored data gives the following relation: log 
(CaT/H/5) (1.33 ± 0.23) log(Fe iiA4570/H;d) - (0.63 ± 0.07). The relation 

between O l A8446/H/5 and Fe llA4570/H/3 is, on the contrary, much more shallow, 
even if the slope is significantly different from 0: log(0 l A8446/H/3) = (0.56 ± 0.14) 
log(Fe IIA4570/H/5) - (0.76 ± 0.05), implying that the systematic O l A8446/H/3 
change is a factor ~ 3.5 for an order-of-magnitude increase in Rreii- 


6.5. Correlation with redshift and luminosity 

Figure [Io| shows that W(0 l A8446) does not correlate with luminosity over a 
range covering 4 dex in luminosity, as defined by our sample and Persson’s data. 
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with the addition of intermediate luminosity sources from Matsuoka et al. (2007). 


The present data do not support earlier claims of a low-ionization Baldwin effect 
in O I, neither as a function of luminosity nor as a function of redshift (i.e., an 
“evolutionary Baldwin effect”). 


The data shown in Figure [T^ are statistically biased in terms of both .2 and 
L (i.e., they are not uniformly distributed over the covered range). In addition, 
there is a less-visible, difficult to quantify “Eddington bias” (see Sulentic et al 


2014 for a discussion on this topic): our high-;^ Pop. B sources may be those with 
the largest L/LEdd possible for Pop. B. Therefore, great care should be exerted 
in interpreting the bottom left panel of Figure [T^ before claiming a correlation, 
even if the correlation coefficient is -0.564, with a probability <10“^ of a chance 
correlation. Accepted at face value, the trend between Mb and W(CaT) would 
imply an increase in prominence of CaT with L and z. We see two possible 
explanations for this result, in addition to an intrinsic luminosity effect: (a) a 
larger ratio CaT/Fe llA4570 in the high-L sample; (b) a systematic difference in 


the Fe llA4570/H/3 distribution between our sample and the sample of Persson 


(1988). 


The correlation between CaT/H/5 and Fe llA4570/H/3 implies that a strong 
CaT emission is associated with a strong Fe llA4570 emission. The bottom left 
panel of Figure shows the i?Feii distribution for our sample and Persson’s. They 
are not significantly different according to generalized Wilcoxon’s tests (repeated 
also without censored data). Therefore, the trend of Figure [l0| cannot be ascribed 
to a systematic difference in i?Feii- The bottom right panel of Figure shows the 
CaT/H/5 distributions. In this case, the same two-sample tests yield a significant 
difference. Figure [11 shows the CaT/Fe llA4570 distributions (only for detected 


data) that are probably least affected by selection effect^ and they are found 
significantly different with a confidence above the 2-cr level: the CaT/Fe llA4570 
ratio average is larger by a factor ~ 1.7 for the ISAAC sample. This effect can 
account at least in part for the L and z trends of Figure and is a potentially 
important result that will be discussed in §7.3| and §7.4[ even if a confirmation 
from a larger sample and a set of observations obtained with the same instrument 


t There is no a-priori reason why we should select sources with a biased distribution of 
CaT/Fe iiA4570 since we are able to measure the ratio CaT/Fe iiA4570 over a broad range of 
values (~ 0.05 — 0.7) and we have undistinguishable itpeii distributions for the Persson and 
ISAAC sample. 
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including low- and high- 2 ; sources is needed. 


Discussion 


7.1. Comparison with photoionization models 


We computed photoionization models using the code Cloudy (Ferland et al 


1998, 2013) in order to constrain the physical conditions of the region where CaT 


and O I A8446 are emitted. The simulations include a 371-level Fe ll ion that allows 
some limited comparison with the observations of optical and IR Fe ll emission. 
Simulations span the density range 7.00 < logun < 14.00 cm“^ and ionization 
parameter range -4.50 < log U < 00.00, in intervals of 0.25 dex assuming plane- 
parallel geometry. The spectral energy distribution of the ionizing continuum is 


the one of Mathews & Ferland (1987) that is considered a standard for quasars. 
Open geometry, line thermal broadening and solar metallicity are assumed. We 


consider column density values 10^^ and 10^® cm“^. Figure 1^ shows isopleths of 
the predicted CaT/0 l A8446 ratio as a function of U and uh for the two values 
of column density. Figure [l3| shows the panels isopleths in the plane (77, uh) for 
CaT/H/5, O I A8446/H/3, CaT/Pa9 and O l A8446/Pa9, again for the two values 
of column density. 

From a first comparison between the photoionization prediction at Z = Zg 
and the observed line ratios we can draw the following preliminary conclusions: 


• The geometrical depth of the fully-ionized zone does not exceed the geo¬ 
metrical depth of the gas slab in the case of = 10^® cm“^. There is a 
depth range where the gas is only partially ionized, even at the highest U 
(log 77 r\j 0). In other words, the ionizing photon flux is not able to make the 
cloud optically thin. In a photoionization context, O l A8446 and especially 
CaT emission will be observed only if the emitting gas remains optically 
thick to the Lyman continuum. Therefore, in the high ionization parameter 
area of the plane (nH,77) we still have some emission of O l A8446 and CaT 
that it is not predicted when W is smaller, 10^^ cm“^. 

• Constraints on the density of the emitting regions follow from maximum 
and minimum measured values of CaT/H/3, O l A8446/H/3, CaT/Pa9 and 
O I A8446/Pa9. Regardless of N^, the ratios involving CaT (when CaT 
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emission is detected at a high significance) suggest that emission occurs only 
at relatively high density log uh > 11.00, with an upper limit to the ionization 


parameter log U 

(poo^. 


-1.5, in agreement with the results of Matsuoka et al. 


• The median values of the CaT/H/5 (~ —0.6) and CaT/0 l (~ 0.2) ratios 
do not favor the area in the plane (17, uh) at very low-ionization and high- 
density (log 12, logl7 < —2.5), unless density becomes log uh ~ 13. In 
this region significant emission of CaT - but not of O l A8446 - can occur. 


• The ratios of CaT and O l A8446 normalized to H/3 and Pa9 look strikingly 
similar (Figure [l3|). H/3 and Pa9 are practically interchangeable. The upper 
half of Table lists the logarithm of the observed flux ratios for H/5 and Pa9 
(BC only), O I A8446 and CaT. 


The observed equivalent width of CaT can be 50 A, and W(0 i A8446) 
< 25 A. These large equivalent widths are observed in both our sample and the 
Persson sample. Reaching values that high for CaT is possible in the framework 
of photoionization within the broad uh lower limit and U upper limit set above. 
The W(CaT) is strongly dependent on the column density in the area centered 
at (nH,C) ~ (11.5, -2.5): the larger column density case yields an almost 10-fold 
increase in W(Ca ll), while the effect is much lower for O l A8446. 


If we want to impose the conditions that O l A8446 and CaT are emitted 
mainly within the same range of physical conditions, then a high metallicity (5 
Zq) seems to be required in order to account for the largest observed W(0 l A8446) 
with a reasonable covering factor (~ 0.2). Metallicity solar or higher than solar 
is a very likely condition for HE sources. They are among the most luminous 
quasars known and there is a well defined correlation between Z and luminosity 
(e.g.. Shin et al. 2013). Among Pop. B sources the HE quasars show large Rpeii- 
Therefore, strong CaT and O l A8446 emitters (the values of Table are loosely 
correlated) would benefit from large column density and high-Z. We remark that 
the condition of high-Z is needed to explain especially the O l strength if we 
require that O l A8446 and CaT are emitted under similar conditions. This might 
not be necessarily the case. 


At variance with CaT, O l A8446 can be emitted over a wide range of den¬ 
sities: 9.00 < logUH A 12.00. We interpret the difference in behavior as due to: 
1) the similarity of neutral Hydrogen and Oxygen ionization potentials (~ 13.6 
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eV) and 2) the Bowen fluorescence mechanism that is strongly influenced by the 
ionizing photon flux. The Bowen mechanism has indeed been found to be the 
major contributor to the O l intensity in most AGN studied by [Rodriguez-Ardila 


et al. (2002a) and Landt et al. (2008). Therefore O l A8446 emission can originate 


in deep regions exposed to a large ionizing photon flux provided that the column 
density is high ~ ]^q 23-25 These properties are consistent with those of 


the VBLR seen in Pop. B sources (Marziani et al. 2010) 


To frame these results on physical conditions in a broader scenario encompass¬ 
ing kinematics, we have to consider inferences that come from analysis of internal 
line shifts and proflle differences. A basic result is that the BLR can be separated 
into two main regions: one emitting mostly LILs and the other emitting HILs. In 
Pop. A sources a blue shift and asymmetry (as in the prototypical HIL C IVA1549) 


indicates that the HILs are emitted within a partially-obscured flow (Sulentic et al. 


2007; Richards et al.|20TT ; Wang et al.|2011 ). Photoionization models further sug¬ 


gest that the HIL emitting region has a relatively low density (uh ~ 10® cm ^), 
column density of cm“^, and a high ionization parameter, U 


10-L 


more 


than 10 times the value appropriate for LIL emission (Netzer 2013). 


The hydrogen lines (typically H/3) have been used as representative lines of 
the LILs. Photoionization models show that LIL emission is associated with a 
density that can be as high as ~ ]^qii- 5 -i 2 column density > 10^^ cm“^ 

and a rather low-ionization parameter logG ~ —2 (Marziani et al. 2010; Netzer 


2013). However, even under these conditions, it is not possible to reproduce the 


intensity of very low-ionization lines, like the Fe II features (Wills et al. 1985). 
Fe II lines need a higher column density ~ 10^^ cm“^) that makes possible an 
extended partially ionized zone at a relatively low electron temperature (T ~ 8000 
K; [Collin-Souffrin et al.||1980] |Collin-Souffrin|[T987t |Joly||1987[ Ferland fc Persson 


1989; Matsuoka et al. 2007). 


In Pop. B the blue shifted component seems to be less prominent. The HILs 


are emitted closer to the central source than the LILs (Peterson & Wandel 1999), 


implying a radial stratification of the emitting region. The distinction BC/VBC 
is the most simple parameterization of the ionization stratification, with the VBC 
being associated with highly ionized (log (7 ~ —0.1), large column density gas 
(log Ac > 23), probably with a broad range of density ([Snedden fc Gask^ 2007 


Marziani et al. 2010)) 
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7.2. A VBC for O l A8446? 


We observe a similarity between the O l A8446 and H/5 profiles (c.f. Landt 


et al. 2008). If we take into account the results of photoionization models it is 


reasonable to investigate the possibility that O l A8446 shows significant VBC 
emission similar to redshifted VBC emission observed in H/5 of Pop. B sources. 


In order to test the presence of a VBC in O l A8446 we added this component 
to the fits for Pop. B sources. The fits for HE0035-2853 and HE1349+0007 are 
shown in Eigure Results are basically the same with or without inclusion of 
a VBC. Even if FWHM O l A8446 and CaT are the same using models with and 
without a VBC, the flux values changed. Inclusion of a VBC reduces the intensity 
in O I A8446 and CaT by 10 - 30% and 20 - 60%, respectively. This reduction 
implies only slightly different physical conditions when the intensity ratios are 
entered in the (V, uh) contour diagrams. The lower half of Table reports values 
for the CaT/H/5, O l A8446/H/3, CaT/Pa9, O l A8446/Pa9, CaT/0 l A8446 ratios 
including a VBC for O l A8446. Comparing the ratios with and without VBC, we 
see that there is no difference within the errors. For example, if we only assume 
a BC, the range of the ratio CaT/0 l A8446 is [-0.59, 0.53] while including the O 
I A8446 VBC implies [-0.62, 0.40]. Therefore, the inclusion of an O l A8446 VBC 
does not significantly modify inferences about density in the CaT and O l A8446 
emitting region: lognn > 11.00 and 9.00 < log nj^ < 12.00, respectively. 

The fits carried out without an O l A8446 VBC indicate that there is no em¬ 
pirical requirement for significant VBC emission to obtain acceptable residuals. 
However the photoionization models indicate that some O l A8446 VBC contri¬ 
bution could be present. Also the shape of the CaT -|- O l A8446 blend bears a 
striking similarity to the H/5 profile in sources where the H/5 VBC is unambigu¬ 
ously present. The S/N of our data are not high enough to constrain any VBC 
emission in xt solutions. We consequently decided to consider a model involving 
a single Gaussian for O l A8446 with eventual consideration for a possible VBC 
contribution. 


7.3. Ca II triplet and O i A8446 (and Fe ii) in the 4DE1 context 

Several important aspects emerge from the previous analysis: in the context of 
a photoionized gas, CaT can be mainly emitted in a low-ionization, dense medium 
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of high column density. In the last decades, there has been considerable progress 
in modelling Fe ll emission in the context of photoionization. As mentioned, the 
same physical conditions suited for Ca ll are also suggested for Fe ll emission. The 
relation between i?Feii and Ca ll/H/5 implies that the same physical mechanisms 
should be operating in the emission of both Fe ll and Ca ll. 

The present study confirms that Ca ll and Fe ll emissions are correlated even 
in very high luminosity objects (Figure [^. Figure shows that our sample and 
the low-luminosity sample of Persson partly overlap and smoothly merge. This is 
another confirmation that the 4DE1 correlations involving low-ionization emission 
line ratios are also orthogonal to luminosity. 


The 4DE1 sequence is, at least in part, a sequence of increasing optical Fe 


II prominence (Marziani et ah 2009, and references therein), and our sample con¬ 
tains mainly Pop. B sources; the Ca ll detection confirms that a dense and low- 
ionization region has to be present also in Pop. B sources. Therefore, the change 
along the 4DE1 sequence should be ascribed to a systematic variation contribution 
of low-ionization region: from a low contribution to Fe II and Ca ll emission in 
extreme Pop. B sources, to a maximum contribution in extreme Pop. A sources 


(Negrete et al. 2012). 


Recent results on the reverberation time lag measured for three Seyfert 1 
galaxies (NGC 4593, Mrk 1511 and Akn 564) strongly support the conclusion that 
Fe II emission originates in photoionized gas located predominantly in the outer 


portion of the BLR (Shapovalova et al. 2012; Barth et al. 2013). Our analysis 


indicates that CaT and O l A8446 can be interpreted as emission from photoion- 


izated gas, with CaT strongly favoured at low-ionization ar 
conditions we suggest are similar to those identified for Fe ll ( 

id high density. The 
Bruhweiler & Verner 

2008; 

Verner et al. 

1999; Sigut & Pradhan 

2003; Sigut et al. 

2004 

). Therefore, 


photoionization models are probably appropriate at least in moderate Fe ll emit¬ 
ters. Refinements in modelling Fe ll emission are needed to reconcile observations 
with theory, in order to reproduce the observed multiplet ratios and specially the 
Fe llf/y/Fe llopt (Sameshima et al. 2011). 


Can also the strongest Fe ll emitters be explained in a pure photoionization 


scheme? According to Collin-Souffrin (1986) and Joly (1987), a different source 


of heating is needed. The strongest Fe ll emitters show continuity with fainter 


Fe II emitters in the 4DE1 sequence (Marziani & Sulentic 2012; Zamfir et al 


2010). By Occam’s razor, invoking an additional mechanism should be avoided if 
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no discontinuity is observed. By virtue of the correlation in Figure strong CaT 
emitter should be also strong Fe ll emitters. If log CaT/H/5 ~ -0.3, we expect 
i?Feii ~ 2. The results of this paper indicate that the strongest observed CaT can 
be accounted for in a photoionization scenario, suggesting that also strong Fe ll 
emitters with i?Feii ^ 2 could be explained by photoionization, if CaT is a valid 
tracer of the Fe ll emitting gas. 


7.4. Implications for Star Formation 


We pointed out (16.5) a possible increase of W(CaT) with luminosity. This 
trend is prone to selection effects and should be viewed with caution. The CaT/Fe 
IIA4570 ratio is, on average, a factor ~ 1.7 larger in our sample than in |Persson| 


(1988). We regard CaT/Fe llA4570 as a reliable measurement. It is unlikely that 
the difference stems from different analysis techniques. In fact our measures are 
corrected for Fe ll and high order Paschen lines, while Persson’s were not. There 
could be several other factors affecting the CaT/Fe llA4570 ratio, for example 


effects associated with microturbulence and spectral energy distribution (Verner 


et al.|2003 ; [Matsuoka et al. |2007| ). The effect of microturbulence on UV Fe ll could 
be especially relevant if Uturb is as high as 100 km s“^. Its effect on the optical Fe ll 
emission could be much less, as tested by CLOUDY simulations for plausible values 
of U and uh (—2.5, 12.00). The possibility of a softer SED has been checked by 


comparing a simulation with the Mathews & Ferland (1987) SED and with the 


NLSyl SED of Marziani & Sulentic (2014). The effect of a softer X-ray spectrum 


is relatively modest, yielding an 25% decrease in CaT/Ee llA4570 ratio at a fixed 
U. This last effect may influence our data since Persson sources are mainly NLSyl 
while our sample includes mainly Pop. B, but is unable to fully account for the 
difference between two samples. 

If 1) CaT/Ee llA4570 traces the calcium abundance as well as the abundance 
ratio of calcium and iron and 2) if calcium abundance scales with other a elements, 
then a CaT/Fe llA4570 increase at higher 2 : could be associated with circumnuclear 
star formation or host galaxy evolution because the ratio [«/Fe] represents a sort 
of “chemical clock” with [a/Fe] higher at earlier cosmic times (e.g., 


Matteucci 


2012, and references therein). An overabundance, relative to Fe, of a elements 


such as calcium is expected if there has been enrichment from a recent burst of 
star formation: the a elements are produced through core collapse supernovae, on 
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timescales < 3-10^ yr, while iron is due to type la supernovae associated with much 
longer timescales, r\j 10® yr dWyse &: Gilmore]|1988t [Matteucci||2003[ ). The sources 
with largest W(CaT) in the Persson| (1988) sample (e.g. the ultra-luminous IR 
galaxy (ULIRG) Mrk 231, Braito et aH2004 ) are also those known to be associated 
with strong star formation and large -Rpeii- These low -2 sources however show 
CaT/Fe llA4570 ratio that is close to the average of the low -2 sample. In the 
simplest scenario of a circumnuclear burst of star formation, enhanced CaT/Fe 
IIA4570 ratio may imply that star formation in intermediate -2 quasars has been 
sustained over a shorter timescale than in the case of the large Rpeii sources in the 
low -2 sample or, more likely, that massive stars are continuously forming. It is 
generally accepted that black hole growth and star formation are closely related 
(e.g. Sani et al. 2010; Kormendy & Ho 2013, for a review). Therefore, the higher 
CaT/Fe llA4570 ratio observed in the ISAAC sample could possibly reflect the 
coevolution of black hole and host galaxy (in the framework of massive black 
hole formation and bulge growth after a wet merger) that is expected to occur at 
intermediate and high 2 and to ultimately power the most luminous quasars. 


7.5. A possible geometry of the BLR 


The narrower line width of Fe ll and Mg ll A2800 indicates that the distance 
from the central continuum source could be larger for these lines than for the 
Balmer lines (Sulentic et al. 2006; Marziani et al. 2013b), if the velocity field 
associated with line broadening is predominantly virial. According to the physical 
condition inference from the CaT, a natural site of emission for LILs is the accretion 


disk that provides a high column density environment (Collin-Souffrin 1987; Collin- 


Souffrin et al. 1988; Ferland & Persson 1989). Physical conditions may not be 


strictly the same for O l and Ca ll emission. O l is more likely emitted in a region 
with physical conditions similar to H/5 (^7.1). 


For example, we may envision a configuration of plane parallel distribution of 
BLR clouds above and below the accretion disk. In this configuration the clouds 
dynamics is dominated by gravity and the emitting gas motion is virialized. Some 
of these clouds might be at the same distance from the ionization source as some 
accretion disk regions, and therefore they may share the same dynamics. 


The Ca ll and O l emitting regions should share a similar ionization status 
since the two ions have similar ionization potentials. This (along with the kinemat- 
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ical similarity) lead to the interpretation that the CaT and O l A8446 are emitted 
from the physically overlapping gas. However, they do not necessarily share ex¬ 
actly the same density and the ionization status. From Figure we see that the 
observed O l A8446/H/5 ratio leads to poor constraints on uh and U. If r is similar 
because of similar kinematics, then there is a large range in uh and U that will 
be consistent with the observed O l A8446/H/5 keeping the ionizing photon flux 
the same for the two lines. This indicates that O l A8446 and CaT may have 
a similar FWHM even if they are emitted in somewhat different conditions. A 
possible model involves an illuminated accretion disk in vertical hydrostatic equi¬ 
librium at each radius: Ca ll may arise with higher eflrciency than O l A8446 from 
the deepest regions of the disk. 


The geometry of the BLR is likely dependent on the source location in the 
4DE1 sequence. The 4DE1 is probably a sequence of Eddington ratio, which de¬ 
creases from Pop. A to Pop. B (Marziani et ah 2001). Non gravitational forces 
such as radiation pressure are expected to be operating at maximum strength 


(Marziani et al. 

2008; 

Netzer & Marziani 

2010; 

Ne- 


grete et al. 2013). Among Pop. B sources, we frequently observe strong redward 


asymmetries in the H/3 profiles, and they are very prominent in the spectra of our 
sources. The physical origin of these features is unclear (e.g Punsly fc Zhang|2010 ), 
and may be associated with a non-Doppler process. If redward asymmetries are 
associated with non-virial gas motions, they may be seen as perturbations whose 
amplitude increases toward the line base of the H/3 profile, and thus toward the 
central continuum source. In the profile interpretation we adopt, we isolate a VBC 
in H/5 profiles that should represent this innermost part of the Hydrogen emitting 
region, the VBLR. The VBLR cannot emit Fe ll or Ca ll because ionization level 
is too high. 


As mentioned, however. Hydrogen and Oxygen have very similar ionization 
potentials but not necessarily the same ionization status, and O l A8446 and H/5 
could be both produced in the inner region exposed to a strong radiation field 
Some O I A8446 emission could be associated with the VBLR. O l A8446 emission 
could be favoured also in an accretion disk plus wind scenario, if the innermost 


regions of line emission are the densest and of highest optical depth (Flohic et al. 


2012). The wind/disk scenario is also consistent with the idea of a mostly virial 


BLR for Pop. B, in which an ionization gradient with distance is responsible for 
differences in line profiles and differences in response times found by reverberation 
mapping (the BC/VBC can be seen as a crude phenomenological model of such 




















ionization gradient). It remains to be seen how much O l A8446 can be produced 
in the VBLR. Theoretical profiles for the different fines will be computed in a 
companion paper in order to ascertain how strong an O I A8446 VBC might be. 


8. Conclusion 

We have analyzed a set of Ca ll IR triplet + O l A8446 observations for a 
sample of high- 2 ; luminous quasars, and reported the detection of CaT emission 
in all sources with exception of two. We carried out multicomponent fits that 
included detailed modelling of all emission features - most notably Fe ll and high 
order Paschen lines - in the range 8000 - 9500 A. The fits allowed us to retrieve 
accurate line widths and fluxes of CaT and O l A8446. The new measurements 
indicate that Ca ll and optical Fe ll emission are likely closely related. A possible 
systematic difference in the CaT/Fe llA4570 ratio (which may also account for an 
apparent correlation of W(CaT) with 2 ; and L) has been found and indicates that 
recent star formation may be driving gas enrichment in the intermediate 2 ;-sample. 

We present a preliminary interpretation of our measurements within a pho¬ 
toionization context. Photoionization models show that Ca ll triplet and O l 
A8446 can be emitted in regions with similar physical conditions, although there 
is a range of ionization parameter and density where only O l A8446 emission is 
possible and CaT emission is negligible. We propose a scenario where the geome¬ 
try of the BLR consists of plane parallel distributions of clouds above and below 
an accretion disk. In this configuration cloud dynamics is dominated by gravity. 
Some of these clouds could be at the same distance from the ionization source as 
high density regions within the accretion disk. The clouds and the emitting regions 
of the disk would share the same dynamics (same FWHM) but differ in physical 
conditions. Much CaT emission could come from the disk, while O l A8446 from 
less dense gas located, however, at the same distance. 

In this paper we have followed a predominantly empirical approach with a 
detailed description of data analysis. In a companion paper we will present addi¬ 
tional data along with a study of the ionization mechanism associated with Ca ll 
triplet and O l A8446 emission. 


The authors thank reviewer for the useful comments to improve this manuscript. 
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Table 1: Basic properties of sources and log of observations 
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Table 3: Measurements of equivalent width, fluxes and FWHM of O l A8446 and 
Ca II triplet 


Object name 

(erg cm“^ A" 


BC 


w 

-b (A) 

F 

(erg cm“^) 

FWHM 
(km s-q 

O I A8446 

HE0005-2355 

0.7 ± 0.03 

7.8 ± 2 (2) 

5.2 ± 1 (1) 

3500 ± 820 (1080) 

HE0035-2853 

1.0 ± 0.02 

9.4 ± 2 (2) 

12.1 ± 2 (2) 

5000 ± 370 (410) 

HE0043-2300 

1.5 ± 0.03 

24.5 ± 3 (3) 

37.1 ± 4 (5) 

4000 ± 360 ^90) 

HE0048-2804 

0.5 ± 0.01 

18.4 ± 1 (2) 

10.1 ± 1 (1) 

4990 ± 270 5oO) 

HE0058-3231 

0.8 ± 0.08 

44.4 ± 8 (8) 

35.0 ± 5 (5) 

4960 ± 100 poo) 

HE0203-4627 

0.8 ± 0.03 

13.6 ± 3 5) 

11.1 ± 3 (3) 

6000 ± 250 ^50) 

HE0248-3628 

1.4 ± 0.02 

12.2 ± 2 5) 

17.3 ± 3 (3) 

3490 ± 260 ^60) 

HE1349+0007 

0.7 ± 0.03 

19.9 ± 4 5) 

14.7 ± 3 5) 

4580 ± 680 (850)* 

HE1409+0101 

1.3 ± 0.05 

15.3 ± 2 5) 

19.3 ± 3 (3) 

3100 ± 310 f(310) 

HE2147-3212 

0.6 ± 0.02 

22.1 ± 6 (11) 

13.6 ± 4 (7) 

4300 ± 860 (1450) 

HE2202-2557 

0.3 ± 0.01 

3.6 ± 1 (3) 

2.0 ± 0.4 (1) 

5810 ± 1060 (1370) 

HE2340-4443 

1.3 ± 0.05 

15.9 ± 1 (2) 

22.0 ± 2 (2) 

3430 ± 220 ^20) 

HE2349-3800 

0.6 ± 0.02 

17.6 ± 4 5) 

11.3 ± 2 (3) 

3480 ± 520 ^20) 

HE2352-4010 

2.6 ± 0.05 

16.6 ± 2 5) 

43.5 ± 4 (4) 

1930 ± no qio) 

Ca II A8498, A8542, A8662 

HE0005-2355 

0.7 ± 0.03 

13.9 ± 8 (11) 

9.1 ± 5 (7) 

4600 ± 1070 (1690) 

HE0035-2853 

1.0 ± 0.02 

27.9 ± 2 (2) 

28.8 ± 2 (2) 

4540 ± 170 q80) 

HE0043-2300 

1.4 ± 0.03 

36.1 ± 3 5) 

53.3 ± 4 (6) 

4000 ± 150 qeo) 

HE0048-2804 

0.5 ± 0.01 

4.9 ± 3 (4) 

2.6 ± 2 (2) 

5170 ± 2400 0200) 

HE0058-3231 

0.8 ± 0.08 

59.6 ± 14 (21) 

45.0 ± 10 (15) 

4910 ± 740 (750) 

HE0203-4627 

0.8 ± 0.03 

33.2 ± 3 5) 

26.3 ± 3 (3) 

5960 ± 530 (650) 

HE0248-3628 

1.3 ± 0.02 

31.3 ± 2 (2) 

42.8 ± 2 (3) 

3990 ± 150 (150) 

HE1349+0007 

0.7 ± 0.03 

24.7 ± 7 (7) 

17.7 ± 5 (5) 

4530 ± 940 (1100) 

HE1409+0101 

1.2 ± 0.05 

19.9 ± 6 (7) 

41.4 ± 8 (8) 

3550 ± 500 (510) 

HE2147-3212 

0.6 ± 0.02 

47.4 ± 4 (5) 

28.1 ± 2 (3) 

3990 ± 150 (150) 

HE2202-2557 

0.3 ± 0.01 

13.0 ± 3 5) 

6.8 ± 1 (1) 

5900 ± 330 (400) 

HE2340-4443 

1.3 ± 0.05 

14.8 ± 10 pi) 

20.1 ± 14 (15) 

3190 ± 1700 (1760) 

HE2349-3800 

0.6 ± 0.02 

19.9 ± 6 (6) 

12.3 ± 3 (3) 

3520 ± 700 (730) 

HE2352-4010 

2.5 ± 0.05 

21.7 ± 3 5) 

55.7 ± 6 (6) 

3080 ± 180 (180) 


Notes. ^In parenthesis, we report the uncertainty taking into account the Fe ii template effect. 
^The equivalent width, flux and FWHM values for O i NC are: I.3±I.O A, 1.0±0.7 erg s“^ cm“^ 
and 550±250 km s“^, respectively. 



Table 4: Measurements of equivalent wi 


Object name 

Fe ^^opt 

w 

(A) 

F 

(erg cm“' 

HE0005-2355 

17.2 zb 2 

51.4 zb 8 

HE0035-2853 

40.7 zb 3 

97.0 zb 8 

HE0043-2300 

20.0 zb 1 

78.7 zb 6 

HE0048-2804 

22.3 zb 2 

15.9 zb 2 

HE0058-3231 

27.7 zb 2 

68.3 zb 19 

HE0203-4627 

14.4 zb 2 

41.5 zb 6 

HE0248-3628 

13.6 zb 0.5 

62.7 zb 3 

HE1349+0007 

18.5 zb 1 

51.6 zb 6 

HE1409+0101 

36.2 zb 2 

203.8 zb7 9 

HE2147-3212 

35.0 zb 4 

85.4 zb 12 

HE2202-2557 

12.7 zb 1 

31.1 zb 3 

HE2340-4443 

15.4 zb 1 

50.4 zb 3 

HE2349-3800 

19.9 zb 1 

35.7 zb 2 

HE2352-4010 

17.6 zb 0.3 

147.6 zb 5 


and flux for optical and NIR Fe ll. 


Fe Unir 


w 

(A) 

F 

(erg cm“^) 

Template used 

20.4 zb 5 

12.1 zb 9 

Theoretical 

34.2 zb 4 

32.5 zb 8 

Theoretical 

2.8 zb 0.3 

3.5 zb 0.2 

Theoretical 

18.8 zb 1 

8.2 zb 1 

Theoretical 

47.3 zb 4 

28.9 zb 7 

Theoretical 

25.1 zb 4 

17.3 zb 4 

Theoretical 

28.1 zb 3 

31.1 zb 1 

Theoretical 

38.5 zb 8 

23.2 zb 6 

Theoretical 

3.4 zb 0.4 

3.6 zb 1 

Theoretical 

11.3 zb 2 

5.6 zb 1 

Semi-empirical 

8.1 zb 2 

3.6 zb 1 

Theoretical 

13.0 zb 1 

15.2 zb 4 

Theoretical 

8.6 zb 1 

4.5 zb 1 

Semi-empirical 

12.8 zb 1 

26.8 zb 3 

Semi-empirical 



Table 5: Observed flux ratios for the BC of H/5, Pa9, O l A8446 and Ca ll triplet 


Object name 

log(CaT/Hd) 

log(0 i/H/3) 

log(CaT/Pa9) 

log(0 i/Pa9) 

log(CaT/0 i) 

No 0 I VBC 

HE0005-2355 

-0.78 

-1.02 

0.11 

-0.13 

0.24 

HE0035-2853 

-0.35 

-0.73 

0.28 

-0.09 

0.38 

HE0043-2300 

-0.67 

-0.82 

0.35 

0.20 

0.16 

HE0048-2804 

-0.99 

-0.40 

-0.39 

0.20 

-0.59 

HE0058-3231 

-0.41 

-0.52 

-0.06 

-0.17 

0.11 

HE0203-4627 

-0.32 

-0.69 

0.86 

0.48 

0.37 

HE0248-3628 

-0.60 

-0.99 

0.46 

0.07 

0.39 

HE1349+0007 

-0.63 

-0.71 

0.45 

0.37 

0.08 

HE1409+0101 

-0.50 

-0.83 

0.49 

0.16 

0.33 

HE2147-3212 

-0.31 

-0.63 

0.97 

0.66 

0.32 

HE2202-2557 

-0.93 

-1.46 

-0.07 

-0.60 

0.53 

HE2340-4443 

-1.05 

-1.01 

0.37 

0.41 

-0.04 

HE2349-3800 

-0.54 

-0.58 

0.49 

0.45 

0.04 

HE2352-4010 

-0.77 

-0.88 

0.35 

0.24 

0.11 

0 I VBC 

HE0005-2355 

-0.94 

-1.10 

-0.05 

-0.21 

0.17 

HE0035-2853 

-0.58 

-0.79 

-0.02 

-0.23 

0.21 

HE0043-2300 

-0.67 

-0.82 

0.35 

0.20 

0.16 

HE0048-2804 

-1.14 

-0.52 

-0.52 

0.10 

-0.62 

HE0058-3231 

-0.74 

-0.52 

-0.40 

-0.17 

-0.22 

HE0203-4627 

-0.70 

-0.75 

0.37 

0.32 

0.05 

HE0248-3628 

-0.60 

-0.99 

0.46 

0.07 

0.39 

HE1349+0007 

-0.74 

-0.89 

0.34 

0.19 

0.15 

HE1409+0101 

-0.67 

-0.96 

0.36 

0.08 

0.29 

HE2147-3212 

-0.45 

-0.71 

0.83 

0.57 

0.26 

HE2202-2557 

-0.96 

-1.36 

-0.01 

-0.41 

0.40 

HE2340-4443 

-1.05 

-1.01 

0.37 

0.41 

-0.04 

HE2349-3800 

-0.70 

-0.65 

0.34 

0.39 

-0.05 

HE2352-4010 

-0.77 

-0.88 

0.35 

0.24 

0.11 
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Fig. 1.— Comparison between Pickles stellar theoretical and observed emission of 
main sequence B (upper panel) and G (middle panel) stars and a giant B (bottom 
panel) star with a black-body spectral energy distribution. The black line is the 
observed emission, the grey line is the theoretical emission and the dashed line is 
the black-body emission. 
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Fig. 2.— Calibrated spectra. Pairs of close panels show the H/3 (left) and Ca ll 
(right) spectral regions before continnnm snbtraction for each object of our sample. 
In both panels the best fit is marked by the red color line and the continuum level 
is marked with the cyan line. Abscissae are rest-frame wavelength in A, ordinates 
are rest-frame specific flux in units of 10“^^ erg s“^ cm“^ A“^. 
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Fig. 2.— Cont 
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Fig. 3.— Quasar spectra after continuum subtraction. Abscissae are rest-frame 
wavelength in A, ordinates are rest-frame specific fiux in units of 10“^^ erg s“^ 
cm“^ A“b The left panels show the H/5 spectral region and right panels show 
the Ca II triplet spectral region. The vertical dashed lines are the rest frame 
H/3, 03, O I A8446 and Ca ll A8498, A8542, A8662. Line coding is as follows. 
Thick black lines: broad components; thin black lines: narrow components; blue 
lines: blueshifted components of H/5; red lines: H/5 and Pa9 VBC; grey fines: Fe 
II contribution; thick dashed black lines: Paschen continuum: thick-brown lines: 
high order Paschen lines; thin-dashed pink lines: best-fit model. Note that the 
Ca II triplet and high order Paschen lines are shown as blends of individual lines 
and not as their sum. 
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Fig. 3.— Cont 
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Fig. 3.— Cont. 
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Fig. 4.— Upper left panel: comparison between theoretical (black line) and semi- 
empirical (grey line) template, as obtained by Garcia-Rissmann et al.j (2012). Up¬ 
per right panel: models of the NIR spectral region for HE1349-I-0007, using the 
two templates. The Fe ll prominent feature at ~ 9200 A is less strong in the semi- 
empirical template, leaving a flux dehcit in the observed blend with Pa9 A9229. 
The lower panels emphasise the difference between the two Fe ll templates (thick 
grey line; left: semi-empirical, right: theoretical). Line coding is the same as for 
the previous figure. 
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Fig. 5.— Stellar contribution of HE2202-2557 using different spheroid masses. 
The thin black and grey lines marked a host galaxy with masses of 10^° and 10^^ 
Mq, respectively. The thick black line marked the host galaxy with a mass of 
10^^ Mq. The grey-dashed line indicate the galaxy with the higher mass 10^^ M©. 
Metallicity 2Zq and age 2.4 Gyr are assumed in all cases. Abscissae are rest-frame 
wavelength in A, ordinates are rest-frame specific flux in units of 10“^® erg s“^ 

cm“^ A“h 
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Fig. 6.— Quasar spectra modeled without Ca ll triplet contribution but will 
all other features present. A positive residual at 8542-8662 A indicates the 
presence of the Ca ll triplet. Vertical dash-dot and dot lines mark O l A8446 and 
Ca II A8498, A8542, A8662 rest-frame wavelengths. The thick black line traces the 
assumed continuum (for HE2202-2557 it includes the best fit host galaxy template) 
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Fig. 7.— Top panels: Histograms of W(0 l) (hatched) and W(Ca ll) (shaded) for 
our new ISAAC sample and Persson (1988). Arrowheads indicate upper limits. 
Middle panels: distributions of W(0 l A8446) (left) and W(Ca ll) (right) in our new 
sample (dotted) and from [Persson (|1988|) (shaded). Bottom panel: distributions of 


Fe IIA4570/H/3 (i^Peii) and CaT/H/3 for our ISAAC sample (hatched) and Persson 


( 1988[ )’s (shaded). 
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Fig. 8.— Top: Comparison between the Ca ll triplet, O l A8446 and H/5 FWHMs, 
in km s“^. Abscissa is the H/3 FWHM (full profile, i. e., BC+VBC), ordinate 
is the Ca ll Triplet and 01 A8446 FWHM. Circles represents the O l data, while 
squares are Ca ll triplet data. Black symbols are from our sample and grey symbols 
represent the Persson sample. The diagonal dashed line is the equality line. The 
vertical-dashed line at 4000 km s“^ indicates the separation between Pop. A and 
B. Middle: same as above but for FWHM H/3bc- Bottom: FWHM(CaT) vs 
FWHM(0 I A8446). 
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log(FeII/HJ) 


Fig. 9.— Top: relation between the observed ratios log(Fe llA4570/H/3) and 
log(CaT/H/5). Circles indicate our sample, while squares indicate Persson sample. 
Gray symbols indicate to Pop A. sources, while black symbols indicate to Pop. B 
sources. The arrows mark upper limits. Trendlines are represented by grey dashed 
lines. Bottom: same for log(0 l A8446/H/5bc) vs log(Fe iiA4570/H/5bc)- 
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Fig. 10.— Equivalent width of CaT and O l A8446 as a function of 2 : (top) and 
absolute B magnitude (bottom). Filled circles: sample presented in this paper; 


|Persson| 

1 

); open squares: data from 

Matsuoka et ah 


(2007). Grey dashed lines are trendlines. 
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□ Persson 


log(CaII/FeII) 


.— Distribution of the CaT/Fe llA4570 for the samples presented in this 
batched) and in Persson (1988) (shaded). 
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Fig. 12.— Isocontours of CLOUDY simulation results for line ratio log(CaT/0 l 
A8446) Abscissa is Hydrogen density in cm“^, ordinate is the ionization parameter, 
both in logarithmic scale. Top: models with Nc of 10^^ cm“^. Bottom: models 
with Nc of 10^^ cm“^. The dashed red line marks the maximum and the minimum 
value from the observations. All our measurements are inside this region. 
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Fig. 13.— Isocontours of CLOUDY simulation results for line ratios log(CaT/H/5), 
log(0 I A8446/H/3), log(CaT/Pa9) and log(0 l A8446/Pa9). Abscissa is Hydrogen 
density in cm“^, ordinate is the ionization parameter, both in logarithmic scale. 
Top: model with column density Nc = 10^^ cm“^. Bottom: model with Nc = 10^® 
cm“^. 
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Fig. 14.— Fits with and without a O l A8446 VBC O l A8446. The upper panels 
show the model with only BC. The middle panels present model fits with O l 
A8446 VBC. Meaning of line coding is the same of Figure In the bottom panels 
the continuous line mark the model with O l A8446 BC only, while the dashed line 
marks the model with BC+VBC. We don’t observe a significant difference. 
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A. Notes on individual sources 


HE0005-2355 - CaT and O l A8446 are completely blended. This is the source 
where Pa7 is detected. Emission on the red side of the Ca ll triplet cannot be 
attributed to high order Paschen lines because the decrease in intensity should be 
exponential and this is not observed. [S ill] A9531 is clearly seen. 

HE0035-2853 - One of the highest S/N spectra. The similarity between H/5 
and the O l A8446 + CaT blend is remarkable. We see significant Fe ll emission 
which is also strong in the optical region. There is a clear asymmetry in the Pa9 
+ Fe II blend at ~9000 A due to Fe ll. We can also see Fe ll on the blue side 
of O I A8446. It is underestimated in both theoretical and semi-empirical Fe ll 
templates. 

HE0043-2300 - The contribution of NIR Fe ll emission is weak (optical Fe ll 
is also weak). The whole emission feature at ~ 9200 A is therefore likely to be 
Pa9. Although the spectrum is rather noisy we can detect high order Paschen lines 
(specially PalO and Pall). Some excess emission on the blue side of O l A8446 
could not be reproduced. 


HEOO 48 - 28 O 4 - This quasar shows the weakest CaT emission, if any. The 
blend fit is dominated by O l A8446. Emission on the blend red side is well fit 
by high order Paschen lines. Although Fe ll is weak in the optical spectrum, it 
is reasonably strong in the NIR. Fe ll contributes ~ 70% of the observed flux in 
the bump at ~ 9200 A. Another interesting feature involves the prominent rise in 
intensity at ~ 8200 A which may be the beginning of the Paschen continuum. An 
I band magnitude value is reported in the USNO-B Catalog (Monet et al. 2003), 
originating from photographic data taken ~ 20 years ago. If we normalize the H/5 
spectrum to the USNO-B I magnitude, we get a significant continuum flux increase 
of 25%, which does not produce a qualitative change in our sample statistics. Due 
to the large period of time between photographic observations and our data, we 
prefer not to re-normalize the H/5 spectrum. 


HE0058-3231 - O l A8446 is clearly detected but the CaT detection is uncertain. 
The higher order Paschen line emission could not be reproduced, especially at ~ 
9000 A. On the blue side of the spectrum the presence of emission that decreases 
towards shorter wavelengths is evident and may hint at the presence of Paschen 
continuum. 


HE0203-4627 - O l A8446 shows a blueshift 




450 km s ^ similar to the 
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one observed in the Hydrogen lines. O l A8446 and CaT are completely blended, 
showing a smooth profile but with a clear red asymmetry revealing the presence 
of CaT emission. The low intensity of Pa9 indicates that the contribution from 
higher order Paschen lines is small. For instance, the contribution from Pal3 to 
the Ca II A8662 profile is <10%. For this reason the emission at ~ 8900 A cannot 
be ascribed to Paschen lines. The Fe ll template underestimates the intensity of 
this feature. Fe ll emission is nonetheless favored because Fe ll is three times 
stronger than Pa9 at 9200 A. 


HE0248-3628 - The high S/N of the spectrum allows us to detect Ca ll A8662, 
but does not allow us to distinguish Ca ll A8498 from A8542. H/3 and NIR profiles 
show a strong similarity. Fe ll emission is twice as strong as Pa9 and PalO and 
can partially account for the observed emission at 8900 A. There is also a small 
feature at ~ 9400 that has been identified with Fe ll ( Garcia-Rissmann et al.|2012 ) 
but that is not modelled by theoretical template. Narrow [S ill] A9531 width is 
comparable to the narrow optical lines. We could not reproduce emission observed 
on the blue side of O l A8446. 


HEI 349 +OOO 7 - We found evidence for a narrow component in O l A8446. For 
many years O l A8446 was thought to be produced only in the broad line region 
(Grandi 1980). More recently, evidence was found that part of this line might 
arise in the narrow line region (Landt et al. 2008). Fe ll (> 50% of flux in the 
Pa9 region) dominates the fit with higher order Paschen lines making a minor 
contribution. 


HEI 4 O 9 +OIOI - The optical spectrum shows prominent Fe ll emission while 
the NIR spectrum may indicate less prominent Fe ll emission. Some Paschen 
continuum emission may be present. 

HE2147-3212 - Low S/N spectrum. The only feature clearly observed is the 
O I A8446 + CaT blend. Ca ll A8662 can be readily seen. We cannot fit the Pa9 
line with good accuracy. 

HE2202-2557 - The only source spectrum requiring subtraction of a host galaxy 
contribution. Host galaxy absorptions are clearly visible (especially the Ca ll 
triplet). The stellar spectrum contributes 50% of the quasar continuum luminos¬ 
ity. After continuum subtraction the O l A8446 -|- CaT emission blend appears 
similar to the H/3 profile. O l A7775 and O l A8446 show the same width, with 
O I A7775 stronger than O l A8446. Emission between these lines can be well 
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fit by the Paschen continuum. Higher order Paschen lines are present but not 
unusually strong. O l A8446 and CaT lines are weak, with equivalent width of a 
few angstroms. The reported flux and equivalent width errors in Table are esti¬ 
mated in a consistent way with the other sources. In the case of HE2202-2557 a 
large uncertainty is associated with continuum placement. Repeating the specf it 
analysis indicates that ± 50% may be a more realistic estimate of the flux and 
equivalent width uncertainty. 


HE 234 O -4443 - Modest S/N at 9200 A prevents a reliable fit of the Pa9 line. 
High order Paschen lines and Fe ll measurements have a large uncertainty. O l 
A8446 and CaT are completely blended and O l A8446 dominates the fit. An I 
band magnitude value is reported in the USNO-B Catalog (Monet et al. 2003). 
The I magnitude would change the flux scale by 40%. As for HE0048-2804 we 
prefer to keep the flux scale obtained from our calibration. 


HE2349-3800 - Low S/N makes it difficult to measure Pa9 and Fe ll. The O l 
A8446 -|- CaT blend is the only distinct feature. We can isolate the contribution 
of the Ca ll A8662 line. There is a bump on the blue side of O l that we cannot 
reproduce. 


HE2352-4010 - In this case we distinguish the contribution of Ca ll A8662 as 
well as the other two lines of the triplet. We observe a bump on the blue side of 
O I A8446 that cannot be modeled. 
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B. Error estimates 

Synthetic Gaussian line profiles were built to estimate statistical errors in¬ 
volved in the line parameter measures. Since we are not dealing with isolated 
emission lines, we have to consider the effect of blending on error estimates. For 
the optical spectra, we considered the broad component (BC) and the very broad 
component (VBC) of H/5. For the NIR spectra, we took into account that O l 
A8446 and CaT are blended together. Model blends were constructed using O l 
A8446 and the three individual components of CaT, to which we added Gaussian 
noise. Models profiles were built considering for the typical FWHM and equivalent 
width measured in our spectra. For H/3 BC we assumed FWHM = 4800 km s“^, 
W = 24, 40, 56, 70 and S/N = 20, 30, 40, 50, 60. For H/3 VBC, we considered 
FWHM = 10500 km s“\ W = 5, 10, 20, 30, 50 and the same S/N values. For 
O I A8446, we considered FWHM = 4500 km s“\ W = 10, 15, 20, 40 and S/N 
= 10, 20, 30, 40, 55, and for CaT FWHM = 4300 km s-\ W = 10, 20, 30, 50, 
60 (sum of the three lines) and the same S/N values. In this way we were able 
to build blended spectra with profiles similar to the ones seen in the real data of 
each quasars. 

Let pq be a blended profile, with a given FWHM and W, to which noise has 
been added. Around it a set of profiles po,ij was built: 1) changing FWHM and 
keeping flux and the peak of the line fixed, and then 2) changing flux (and W) but 
keeping FWHM fixed. yAj values were computed for fixed noise level as a function 
of FWHM and flux. The original profile Pq is associated with the minimum 
Xu,min- Errors are set through the ratio Fy = precisely by values 

of Fy that correspond to the probability 0.32 (1-cr confidence level) of reaching 
a yAj value because of statistical errors (i.e., associated with noise). Values of Fy 
are computed for the degrees of freedom set by the number pixels in wavelength 
covered by the lines, minus the number of parameters appearing in the line model. 
The Fy values provide confidence intervals as a function of FWHM, W, and S/N 
that define our uncertainty range. The FWHM error was taken from models where 
the FWHM varies and the the peak of the line is fixed. Flux can be written as 
the product FWHM times peak intensity. Error propagation on FWHM and peak 
intensity yields the total error on flux. The left panel of Figure [T] shows an example 
of blended synthetic profiles with FWHM(0 l) = 4500 km s“^, FWHM(Ca ll) = 
4300 km s“^, W(0 l) = 20 A, W(CaT) = 20 and S/N = 20. The right panel 
shows the uncertainty range on FWHM and flux estimated from the appropriate 
Fij value. The parameter F = I Xu,min shown in Figure as a function of 
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FWHM and intensity of CaT and O l A8446. F ~ 1.03 provides 1-a confidence 
ranges on the parameters. This is basically the classical technique for error 
estimates described by Bevington & Robinson (2003) and by Press et ah (2002). 


We tested the effect of Ca ll+O l blending also with Monte Carlo simulations 
that allowed us to simultaneously fit flux and FWHM. The cases studied for the 
xl analysis were all repeated for different values of S/N ratio. The Monte Carlo 
simulations involved the repetition of the blend fits for a large number of trials 
(!5i 1500 for each case) with different realizations of Gaussian noise patterns. The 
uncertainty on line parameters was then retrieved by the distribution of the mea¬ 
surements obtained changing the random noise patterns. We found that the errors 
are correlated with the ones estimated through the xl analysis, but are signifi¬ 
cantly smaller, between a factor 1.3 - 1.5 and a factor > 2 for the flux of CaT. 
The lower errors can be understood by considering the nature of the Monte Carlo 
Method and by the fact that we are dealing with spectra whose S/N is relatively 
high, > 20 save in one case. In these conditions, the fitting routine converges to 
values close to the “true” ones i.e., the starting values in absence of noise, yielding 
a rather peaked parameter distribution for different noise realizations. 


In addition to Gaussian noise, another relevant source of error is continuum 
placement. To estimate its effect on ffux, we considered a continuum intensity 
displacement of ±lcr with respect to the best fit. The W uncertainty was then 
estimated by propagating the errors for continuum and flux following the definition 
of W. 


The uncertainty due to the inclusion of the Fe ll template has been esti¬ 
mated through Monte Carlo simulations considering all components assumed to 
be present in the whole spectral range 8000 - 9900 A. The flux and width assumed 
for each fine are the ones reported in Table 2 and 3 i.e., the ones measured on the 
observed spectra. We run two cases for each object, in which the fitting routine 
(1) considers only one template as listed in Table 4, and (2) is left free to choose 
between the two different templates (theoretical or semiempirical). The ratio be¬ 
tween the errors found for case 2 and case 1 gives us the uncertainty added by 
Fe II over the uncertainty associated with other components. We multiplied xt 
uncertainties by this ratio (always > 1) to include the effect associated with the 
Fe II template. The uncertainty increase ranges from a median of ~ 10 % to a 
maximum 50%. The effect depends on S/N and W and, in the Monte Carlo 
simulation context, should be ascribed to the inability of the fitting routine to 
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identify which template is the best. The uncertainties computed with and without 
the effect of the Fe ll template are reported in Tables and Figures and 
error bars include Fe ll template uncertainty. 
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A (A) 


Fig. 1.— In the left panels we show the synthetic prohles for the blend O l A8446 + 
CaT with FWHM(0 l A8446) = 4500 km s-^ FWHM(Ca ll) = 4300 km s-\ W(0 
I A8446) = 20 A, W(Ca ll) = 20 and S/N = 20, 40 (top and bottom, respectively). 
The dark line marks the model with the minimuum while the dashed line in 
colour indicates the model at 1-a for FWHM and Flux (yellow and orange). The 
vertical dashed line are the rest-frame for O l A8446 and Ca ll triplet. In the 
right panels the behavior of the F statistics is shown for the 2 different S/N values 
(pale: S/N = 20, dark S/N = 40). Contour lines (traced on the 3D paraboloids 
and on the xy plane) trace the 1-a confidence ranges on FWHM and flux reached 
for F 1.03. 




































